EPJ manuscript No. 

(will be inserted by the editor) 



g Mechano-chemical coupling in growth process of actin gels and 

o 

(N 

^ a symmetry breaking instability 

, Ken Sekimoto"'^'*, Jacques Prost^'', Frank Jiilicher^'^, Hakim BoukcUal^ and Anne Bernheim-Grosswasser^'* 
Physico-Chimie, UMR168 Institut Curie, 26, rue d'Ulm 75248 Paris Cedex 05, France 

Received: date / Revised version: date 



C/3 

O 



> 

l> 
O 

o 
cn 
o 

O 
• 1— I 



X 



Abstract. It has been observed experimentally that the actin gel grown from spherical beads coated with 
polymerization enzymes spontaneously breaks the symmetry of its spherical shape, and yields a "comet" 
pushing the bead forward. We propose a mechano-chemical coupling mechanism for the initialization of 
this symmetry breaking. Key assumptions are that the dissociation of the gel takes place mostly in the 
region of the external surface, and that the rates of the dissociation depends on the tensile stress in the 
gel. We analyze a simplified two-dimensional model with a circular substrate. Our analysis shows that 
the symmetric steady state is always unstable against the inhomogeneous modulation of the thickness of 
the gel layer, for any radius of the circular substrate. We argue that this model represents the essential 
feature of the three-dimensional systems for a certain range of characteristic lengths of the modulation. 
The characteristic time of the symmetry breaking process in our model depends linearly on the radius 
of curvature of the substrate surface, which is consistent with experimental results, using spherical latex 
beads as substrate. Our analysis of the symmetry breaking phenomenon demonstrates aspects of mechano- 
chemical couplings that should be working in vivo as well as in vitro. 

PACS. 87.17. Jj Cell locomotion; chemotaxis and related directed motion - 87.15.Rn Reactions and kinet- 
ics; polymerization - 62.40.-|-j Anelasticity, internal friction, stress relaxation, and mechanical resonances 
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1 Introduction 

Polymerization of actin is one of the main mechanisms 
responsible for cellular motility. Filaments of F-actin are 
polymerized on the cytoplasmic side of a cellular mem- 
brane with the barbed ends oriented towards the surface 
of the membrane. The branching of the actin filaments 
takes place mainly in the vicinity of the surface. The re- 
sulting branched F-actin filaments take a form of a soft 
elastic solid, Tl, which we call an actin gel or a gel, sim- 
ply. This actin gel pushes the cellular membrane outwards. 
Polymerization of actin gels is also a locomotive mecha- 
nism for intracellular bacteria like Listeria monocytogenes^ 
and perhaps also for the endosomes and lysosomes ( 2 ). 
In these cases, the gel is grown in the form of a comet. 
This comet pushes the bacterium forward. For the pur- 
pose of understanding this mechanism of motility, various 
experimental model systems have been developed using 
both biochemical and biophysical approaches. 

Biochemical approaches have isolated the basic cyto- 
plasmic ingredients needed for the motility of Listeria 
monocytogenes [HlKl) E^ctin and ATP for the formation of 
F-actin filaments, (2) Arp2/3 as the cross-linker and/or 
the nucleator of the F-actin growth (i.e., the precise role 
is still under debate), (3) ADF as the depolymerization 
factor at the pointed end of F-actin, (4) the capping pro- 
tein, (5) a bacterial protein called Act A expressed on the 
surface of Listeria which is necessary for inducing poly- 
merization from the surface. These ingredients constitute 
a model cytoplasm for the motility. 



il coupling in growing actin gels 

Biophysical approaches have taken the Listeria as a 
model system of cellular motility. Furthermore, a bio- mimetic 
in vitro system of the bacterial motility has been intro- 
duced. This system consists of a spherical latex bead, 
coated by the enzymatic protein complexes, Act A ( 0], 
or a fragment of its homologue from human cells, 
WASP [Wiskott-Aldrich-Syndrome Protein] [HIE]. The cy- 
toplasm has also been replaced by the reconstituted cyto- 
plasm 7 . Despite the spherical form of the bead, the gel 
has grown in a shape of a comet, like the bacteria Liste- 
ria (see for example Fig. 2 of [5 ). Fig^shows the initial 
stage of the creation of the comet, observed using fiuores- 
cent probe attached to actin monomers. 
This process bares the signature of a spontaneous symme- 
try breaking, which is the subject of the present paper. 
The phenomenon of the symmetry breaking is relevant to 
some biological systems of sub-cellular level. For instance 
the endosomes, which consist of spherical soft substrate 
(liquid vesicle), grow a comet Also, the motility of 
a mutant Listeria, which moves preferentially in lateral 
directions 8 grows the actin gel by breaking its cylin- 
drical symmetry. Our principal aim is to assess, through 
the study of the symmetry breaking, the relevance of the 
elastic aspects to the biological motilities based on the 
polymerization of protein filaments, and to provide for 
several basic ingredients related to the mechano-chemical 
coupling. The three ingredients essential for explaining the 
symmetry breaking are (details will be given in sections 
§2-§4): 
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Fig. 1. Experimentally observed time sequence of the actin gel g 
is 10 fim. The bead is coated with a fragment of WASP, and it is 
actin monomers in the gel is visible by a fluorescent marker. The 
gel has started to grow. 

(i) The creation of a tensile stress due to the curvature 
of the substrate surface {^i'- 

As the gel is continuously created at the bead surface 
(at the radius, r — tq), the part that has been aheady 
formed is continuously pushed outwards (r > tq). Since 
the perimeter (27rr) increases as r, and since the surface 
has a closed topology, the gel is stretched by the ratio, 
r/ro(> 1). 

(ii) The concentration of the tensile stress by a geo- 
metrical effect O: 

The gel layer around a bead is in mechanical equilibrium, 
so that the integrated tension across the layer thickness of 
the gel must be constant along the surface of the bead. In 
particular, if the thickness is locally thinner, such a region 
must bear a stronger tension in order to support the same 
integrated tension. 

(Hi) The acceleration of dissociation of the gel under 
tensile stress (2J: 

We suppose that, under tensile stress, the gel dissocia- 



;rown around a spherical latex bead. The diameter of the bead 
placed in a reconstituted solvent as described in the text. The 
observation started at t = after several tens of seconds when 

tion is accelerated through the mechano-chemical cou- 
pling. This dissociation may be either through the un- 
binding of the branching points along actin filaments, or 
through the depolymerization of actin filaments. 

These three ingredients (i)-(iii) constitute a positive 
feedback loop leading to an instability of the symmetric 
shape of the growing gel. This will be described in ^5.11 
In short, the region of gel with smaller thickness becomes 
preferentially dissociated due to the higher tensile stress, 
implying further thinning of that region. 

All symmetry breaking models take the mechano- 

chemical coupling into account. Previous models have fo- 
cused their attention on the compressive force acting on 
the actin filaments at the polymerization sites, that is on 
the substrate surface. Our analysis takes into account the 
global stress distribution. Of particular importance is the 
tensile stress generated at the outer surface of the gel. In- 
deed, on general grounds, the depolymerization rate must 
be an increasing function of the tensile stress. We show in 
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the following that it leads inevitably to symmetry break- 
ing. In the discussion section ^ after a brief summary, 
we compare, in more details our analysis with the existing 
ones, and suggest experiments designed for distinguishing 
between the different possibilities. 



2 Distribution of stress within the gel with Fig. 2. The cross section of F-actin gel around a bead. The 

relevant stress components in the gel are schematically shown. 



symmetric shapes 



The gel occupies the space between the radii r — rg and 
r — ro + h. The compressive component of the stress at the 
substrate surface (r = ro), (Jrr\^^^, and the tensile component 
at the outer surface (r = ro + h), (t_l_l|^^^^, are indicated by 
the pairs of oppositely oriented open arrows. In this symmetric 
state, the tensile component at the substrate surface, (t±±|^^, 
as well as the normal compressive component at the outer sur- 



face, o-rrl^jj+ft, vanish. 



Suppose that a gel has been polymerized steadily from a 
substrate surface of either spherical or cylindrical shape 
with radius ro, until the gel forms a layer of a thickness h, 
enclosing the substrate surface and keeping its symmetry 
(see, Fig|2l). As already noted, the part of the gel that 
has been formed has been continuously pushed outwards. 
An element of the gel at the radius r is then stretched by 
r/ro times relative to the native state of polymerization. We remark that the present approximation ignores the 
To know the tensile stress in the lateral direction, <j^^ , ""adial deformation due to the lateral stretching, in other 
let us use the "stacked rubber band model" A ^ vanishing Poisson ratfo. Although this 

freshly cross-linked gel at the latex surface (r = ro) is ^^^^ justification for actin gels, the main results of the 
unstretched and has no lateral stresses, <j^^\ = 0. As P^^e^^^^^ P^P^r do not depend on this property. More re- 
the layer is pushed outward, its circumference increases ^^^^^ calculation confirm the validity of this statement 
which introduces a lateral stress (JI]), 



<7±±l = B^-^, (1) 



The shear component of the stress ar± vanishes every- 
where, for symmetry reasons: 



with B being the Young modulus. In particular, when the ar±_ = 

thickness of the gel layer is h, the tangential stress at the 
outer surface of the gel is 



Radial force balance requires that the radial stress, ar 
at the radius r obeys the following equations: 



a^±L+H = B^. (2) L^(r'a„)^-a^^ = (3) 
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for a spherical surface ([HI, see Appendix A), or corresponds to the constant function, h{uj) = h*, with a 

9 constant thickness h* . 

— [rCTrr) - <T±± = (4) 



for a cyhndrical surface [Ilj . Since no external force is 
applied on the outer surface of the gel layer, the normal 
stress must vanish: 



The analysis of the thickness perturbation is done in 
the following two steps: In the first step, we suppose that 
this function, is slightly perturbed from a constant 

function, but we still do not allow for the displacement of 
arr\rg+h = 0. (5) the gel. In the second step, we let the gel layer relax until 

it reestablishes the mechanical balance. We calculate how 
the stress in the gel is distributed in this new balanced 
state. To avoid any confusion, we stress that, the pertur- 
bations {h{Cj) — h* 7^ 0) at the end of the first step does not 
imply the swelling or deswelling of the gel layer. The per- 



Under this condition, the normal stress at the substrate 
surface, CTrrl^Q, can be calculated in terms of the lateral 
stress, a^^_. In the case of cylindrical substrate, we inte- 
grate Eq.Q) from r = rg to r = ro + /i, and have 

T 



(6) 

ro ' turbation rather implies that there is more or less material 

where T is the integrated tension across the symmetric gel of gel along the direction uj than the average. (It could be 

slab, defined by due to the enhanced/depressed polymerization, or, to the 



T= / a^^dr. (7) 

J ro 

Using Eq.([3) we find T = Bh^ / (2ro), and thus (Trr\ro = 
—Bh^/{2rQ). For the spherical substrate, the relation is 
not as simple as the cylindrical case. Still, arr is given as 
an integration of 



depressed/enhanced dissociation of the gel along this di- 
rection.) As we discuss a situation such that there is no 
external force on the outer surface, we require stress-free 
conditions on the outer surface of the gel layer. 

= 0. (8) 



What we demonstrate is that, under the above condi- 
tions, the tensile stress a±±_ mider the reestablished me- 
chanical balance is locally enhanced in the thinned region 
of the gel layer, that is, in the zone of the orientation uj 
In this section, we consider how small perturbations to that satisfies h{uj) < h*. On the one hand, the physical 



3 Concentration of the tensile stress under a 
modulated surface profile 



the surface profile of the gel layer lead to the redistribu- origin of the stress concentration is quite simple and uni- 

tion of the stress components within the gel layer. We versal. In fact the authors have noticed, after completion 

introduce the function representing the thickness of the of the present work, that essentially the same mechanism 

gel layer, h{uj), with the variable lu representing the orien- of stress concentration has been discussed long before in 

tation from the origin. A spherically symmetric gel layer the context of crystal growth under stress (see, for exam- 
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pie, a concise review on the related history in the litera- 
ture ^5). In Appendix B we describe the basic mecha- 
nism of this phenomenon by using an illustrating example 
in a very simple geometry. On the other hand, the direct 
analysis of the present case with the distributed thickness 
h{uj) is difficult, because of the three spatial dimensional- 
ities and the tensorial character of the stress associated to 
this space. We can avoid, however, this difficulty by the 
following lines of reasoning. 

1. We limit our concern to the modulations — h* 
whose characteristic wavelengths are comparable to the 
average thickness, h*. Experimentally, h*/rQ is at most 
about 0.2 jJS]. The radius of curvature of the outer sur- 
face (~ ro + h*) is, therefore, not appreciable in view of 
such short wavelength of modulation. We may then ignore 
the effect of a specific curved geometry of the substrate 
surface except for the fact that the curvature gives rise to 
the lateral tension (Tj_j_ in the gel layer. 

2. We notice the following fact: As far as the stress dis- 
tribution inside the gel layer is concerned, the influence of 
the surface profile perturbation is practically limited to a 
region near the outer surface (see, FiglSj. More precisely, 
if the perturbation is characterized by a wavelength. A, 
then the thickness of the disturbed region is also charac- 
terized by A. (The boundary condition far from this layer 
is therefore irrelevant to this disturbance.) For detailed 
arguments, see Appendix C. 



Fig. 3. Schematic representation of the outer gel surface before 
(a) and after (b) perturbation. A perturbation of the surface 
profile with characteristic length A affects the stress profile only 
within a "skin layer" of thickness A. 



spherical one, with (b) a, "slip" boundary condition on the 
substrate surface, to see how the stress in the actual three- 
dimensional case is distributed after the reestablishment 
of the mechanical balance. Aloreover, (c) the neglect of 
the shear stress components within the gel layer is justi- 
fiable for the experimentally realized situation where the 
mean thickness of the gel layer h* is much smaller than 
the radius tq. Wc will formulate these assumptions in more 
details below: 

(a) We consider the gel layer grown around a two dimen- 
sional circle of radius tq . 

In two dimension, we represent the thickness of the gel 
layer by h{9o) as a function of the angle 9q with < 
9o < 27r, instead of h{6j) above. (See, FigHl h{6o) is 
defined before the reestablishment of the mechanical bal- 
ance.) The lateral components of the stress, which we have 
3. With our wavelength choice in 1., and with the fact denoted symbolically by _L, corresponds now to the az- 
just mentioned above 2., we can justify the study of (a) imuthal direction. We then use the suffix 9 in place of 
a two dimensional circular geometry rather than the real _L hereafter. For example, we write are for ar_L, and agg 
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(Remark: Note that we do not claim this boundary con- 
dition to be always realistic. We rather use this condition 
since it is justifiable for the calculation of the stress distri- 
bution under the modes of perturbations with q ~ 
See the argument 2 above and the appendix C for the de- 
tails. ) 

(c) The shear stress <7re within the gel layer is negligible. 

Fig. 4. Definition of the angler variable 9o and the hight func- 

As mentioned above, the experimental value of h*/ro is 

tion h{9o) of a gel with modulated thickness due to depolymer- 

<C 1. In such situation we may, in the lowest order ap- 

ization before an elastic deformation reestablishes mechanical 

proximation, estimate the magnitude of the shear stress, 

equilibrium. 

with a parabolic profile of the shear stress are- are — 
instead of a^^. For smaU perturbations of the thickness, /ie^(ro + h - r){r - ro)/r^ for r^ < r < r^ + h, which 
1/1(6*0) - h*\/h* < 1, we may use the linear analysis. Then satisfies the boundary conditions, Eqs.© and Here, 
it suffices to consider the form /i is a constant proportional to the shear modulus ^ of 



h{eo) = h*[l + eqcos{qeo)], (9) 



the gel. The magnitude of (7^9 is, therefore, at most of the 



order of jleq{h* / r^f . We compare this with the change of 
where the integer q indicates the number of nodes of the perturbations of the thickness, which is of 

spatial undulations, and is supposed to be small (|e,| « ^qBh*/ro (see Eq.©). Then <Jre is smaller than this 

1). The characteristic wavelength for the the 9-th mode is a factor of h*/ro, and is therefore negligible, 
about 2^ro/g, and the restriction (1) is represented as ^^^^ Appendix D we show how the stress distribu- 

^ 27r7-o tion within the gel is calculated for the model described 

^ — h* ■ 

(Remark: Besides our purpose of analysis, the two-dimensiouifey ^^1°^ ^^^^ o^l^ the resuhs for the ten- 

geometry applies rather directly to a Listeria mutant H ^^^^ ^t^^^^^s ^sel.+h at the external gel surface and the 
mentioned in 32 This mutant moves preferentially in lat- "^o™^! compression, - <Jrrl, at the substrate surface: 

o 2 

eral directions, breaking its cylindrical symmetry.) cggj^,^^^ — (cr6ielro+;i)g ^0 ^ ^ 9 , . -e,cos(9^o) + 0(e,') 

(b) On the substrate surface (r = ro), the shear stress is (11) 
negUgtble. -a^^l^-^-arr\rJ,^^, = 0{e,'), (12) 

The slip boundary condition for the shear stress is written ^.^^^ ^ ^ r^^^ bracketed terms with the subscript 

^ e„ = are those terms without perturbation: 



are\r,=0- (10) ('^eelro+Oe„=o = ^^a'^d(~a„|^J^^^p = fx2(seeEqs.@ 
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and In Ea. p2|l . 0{eq^) indicates the terms of at least assumption: (i) On the substrate surface, the polymeriza- 
second order of eg. Since cos(g 0o) = ih{Oo) — h*)/h*, the tion is the dominant process. The pre-exponential factor 
minus sign on the right hand side of Eq. implies that kp represents the kinetic constants in the absence of com- 
the lateral tension is augmented, (Jeelrg > i'^eelro) ^ -o' ™ pressive stress {(Jrrlra ~ p depends on the concen- 
thc thinned portion of the layer, h{9o) < h* . tration of actin monomers in the solvent. In our analysis 

we assume this to be constant. The exponential factor 
4 Mechano-chemical coupling: Growth and represents the fact that the polymerization is decelerated 



dissociation of gel under stress 



In this section, we consider the time evolution of the thick- 



by the compression, o'rrlj.g (< 0)- The parameter Cp has 
been introduced so that —Cp cr,.r|^^ accounts for the in- 
crease in the polymerization potential barrier due to the 

ness of the gel layer. We denote the profile of the thickness , • , u i ^ i, -^i, i 

s> J t' QQgf^ elastic energy (divided by kbJ ) to push out the gel 

at the time t as h{0o, t). We are interested in the chemi- ^^y^^, ^^^^^^^ ^^^^^^ compressive stress. We have ne- 

cal processes which take place on time scales much larger ^^^^^^^ ^j^^ dissociation of the gel at the substrate surface, 

than the establishment of the mechanical balance within ^^^^ ^^^^^^^ ^^^^^ ^^^-^^ incorporated in the model 
the gel. We suppose that the relevant microscopic chemi- j.^^j^ ^^^^^ ^^^^^^^ ^j^^ 

cal processes are the polymerization and branching of the j^^^^^ polymerization of microtubules, it has been 

actin filaments to form the gel, and the unbinding of the ^^^^^ ^^^^ negligence of the depolymerization on the 

branching points and/or through the depolymerization of ^^^^^j^^g ^^^^ pl^^ ^^^^ ^ ^^^^ approximation ^7\. 
actin filaments to dissociate the gel. We adopt a simplified 



version of the model proposed previously |15ll5] : 



dhi9o,t) 

— a 



kpe'^'-^-^-o"" - fc^e'"'"'l'-o+M«o.t)'='' 



(13) 



The second term in the square bracket on the r.h.s. 
of Eq^l represents the gel dissociation. Wc have intro- 
duced the assumption: (ii) The dissociation process is al- 



dt 

where a, hp, kfi, Cp and are positive constants. The pref- most localized on the outer surface of the gel at r = r^ + 

actor a outside the square bracket on the right hand side h{do,t). The pre-exponential factor, —akd, therefore rep- 

(r.h.s.) is a length of about the size of an actin monomer, resents the rate of thickness decrease which occurs due 

This represents the rate of conversion between the chemi- to the dissociation of the gel under the stress-free condi- 

cal processes and the change of the thickness, ft. The other tion, '^90\rg+h{eo t) ~ ^- (Remark: Here we can identify 

parameters are described below. o'e6i|r(,+?i(e(, t) ^he tensile stress along the tangent of the 

In the square bracket on the right hand side (r.h.s.) of outer surface, since the correction is of second order of the 

Ea ll3l the first term represents the polymerization at the deviation angle, (ro + h), between the tangential di- 

substrate surface (r = ro). Here we have introduced the rection and the azimuthal direction.) The exponential fac- 
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tor of this term represents the fact that the dissociation physiological concentrations of the actin monomers and 

is accelerated by the lateral tensile stress cree\,.g+h(eo t) "-"^ ^^'^ cross-linker molecules, diffusion does not limit the 

(> 0). The parameter Cd has been introduced so that thickness evolution, h(9Q,t). 

Cd cree\rg+h(eo t) accounts for decrease in the depolymer- The evolution equation Ea. H13|) has a solution corre- 

ization potential barrier due to the release of the elastic sponding to the symmetric stationary state, h{9o,t) = h* 
energy (divided by fceT) when the gel is dissociated un- If we restrict our analysis to the circularly symmet- 

der the tensile stress. We have neglected the dissociation ric profiles, h{dQ,t) = h{t), this solution is stable. In fact, 

of the gel occurring inside the gel. There are good reasons substituting the form h{6Q, t) = h* into Ea. (|13|l . we obtain 

to belive that the gel dissociation is strongly accelerated the equation for x = h* /tq as 
under tensile stress QHI j as compared with spontaneous Cp 9 2 f k„\ 

dissociation under the stress-free condition. In fact, the 

r 11 11 , , ,1 1 ,1 This equation has a positive, therefore physically meaning- 
experiments usmg the full cell extract as the solvent have j- : j- j j o 

, ,1 , ,1 ,, • , r ,1 11 1 ful, solution for k„/kd > 1, Furthermore, if we substitute 

shown that the mean thickness ot the gel layer around f/ > 

the latex bead is much smaller than the average length of ^^"^ form 

the comet produced by Listeria of similar size. It implies h{6Q,t) h [1 + eo(i)] i (15) 

that the Boltzmann factor of the form, e'''''^'"', is crucial to the Eq.l^JJ reduces, up to linear order in eo(t), to the 

determine the dissociation rate. As uee is largest on the following equation: 
outer surface of the gel, we suppose that the gel dissocia- deo{t) eo{t) 



dt tq 

tion occurs mostly in the vicinity of the outer surface. 

with To = kd^^f2o^^ro/a, 

The kinetic equation Ea. (|13() also assumes the follow- 



(16) 



f2o = CdBe--^^ ( 1 + -X ) • (17) 



ing: (Hi) The diffusion of actin monomer is fast enough. \ Cd 

This limits our analysis to a bead radius range smaller Equation l|16|l shows, as already mentioned, that the steady 

than a cross-over size, Tc, separating a stress governed state solution h{9o,t) = h* is stable with respect to per- 

regime from a diffusion controlled regime. Indeed, on the turbations keeping the overall symmetry. This result is un- 

substrate surface, the actin gel is formed from the adjunc- derstandable since a radius with h > h* {< h*) would lead 

tion of actin monomer molecules. And for these molecules to an increase (decrease) of both (— cTrrlrg) and cr^_L|^^_|_,j, 

to reach the substrate surface, they have to diffuse through and these in turn make the r.h.s of Ea. (|13|l negative (pos- 

the network of the actin gel. Previous experimental and itive), leading to a decrease (increase) of h toward the 

theoretical analysis indicates that, as far as the diam- stationary value h*. From Eas. H14|l and (|17(l . x and tq are 

eter of the latex bead is less than about 5 /im, and under functions of three parameters, CdB, — , and Note that 
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To is a few orders of magnitude larger than the microscopic our analysis predicts that a quasi symmetric steady state 
time kd~^ , with ro/a being of order 10"^ and of order should be reached significantly earlier than the onset of 
10. symmetry breaking. This is indeed what is observed. 

(ii) The characteristic time of the instability is propor- 
5 Result tional to the radius of the substrate if the other parameters 

are fixed.: It is reasonable to suppose that Tq represents 

5.1 Symmetry breaking instability 

the characteristic time of the growth of the perturbation. 
We now consider symmetry breaking perturbations and Then, from JTHIl, is written in a scaling form: 

we assume the following form for the gel layer profile: Tq _ ^ -i n 

kd a 

h{6o,t) — h [1 + eq{t) cos(g6'o)] , (18) -^^rhere k^^^ and a play the role of intrinsic timescale and 

with q^O. Substituting the expressions of the stress com- lengthscale, respectively. The dimensionless constant of 
ponents, Eqs.CH and CEJ into Eq.^, where is re- proportionality, nq-\ depends on the properties of the 
placed by we have the following equation up to the ^nd of the solvent through the parameters, c^B, ^, 

linear order of e (t) fe' ^^^^ that, in fact, the quantity Tg thus defined 

shows no dependence on 0), as f2q does not. This 



de.jt) ^ eqjt) 

dt Ta 



(19) 



apparently anomalous behaviour should not be taken se- 



with Tq = tqIj^, riously, because the range of wavenumber validity of of 

„ Y our analysis is limited to o ~ ^r?-.) 

2 + X 

Quantitatively, we can evaluate the characteristic time 
Tq using the experimentally known data in the literature: 

The stationary velocity Vgci at which the gel material moves 

Note that since /2oJ7o ~ 2/y ~ 10, tq ^ t„. Eas.ltT^ and ^ i • • i ^-i^ i r rm 

^ ' y 1 ■< t outward is identified from Eq. as 

(1201) imply the following characteristics of the symmetric , , r i \ c^aaa] , , / ; \ ^.r^ ^ 

Wgci/ro — kd{a/ra)e ""i>-o+h(e(,.t) — kd{a/roje . Com- 

stationary state ft.(0o, i) = /i*. ■ ^^ ■ ^^ ■ r • ir am 

paring this with the expression of Tq obtained from Eas. H20|l 

and lEIl), = kd{a/r^) e""^^ CdBx/{2 + x), we see that 



where x — h* /ro as before. (Remember that x can be ex- 
pressed in terms of the parameters CdB, and ^.) 



(i) The symmetric stationary state is unstable against per- 
turbations which break the symmetry, since all Tq are pos- 



""2 + ^ (22) 



itive. In fact, the applicability of our model is guaranteed '^'goi CdBx 

only in the range of q satisfying q ~ ^lla (gee, . Nev- As described in ^ , CdBx is the decrease in energy barrier 

ertheless, the presence of an unstable mode is sufficient (in units of k^T) in the dissociation of an actin filament 

for the proof of instability. Note also that since tq ^ r^, under tensile stress, compared to the unstressed case. For 
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the effects described in this manuscript to be observable, small but finite source term, e^, 

this decrease must be of order one. Noting that x 2, ri (+\ ^ 

^ = 1M0+6^). (23) 

our analysis requires the combination TgUgoi/(2ro) to be ? 

of order one. The experiment gives Tgym/ro ~ 5 min//im We can solve this equation with the initial condition (0) = 

and Vgci — 1/im/min (note that it is the polymerization 0: 

rate under stress), which leads to TgWgei/(2ro) ~ 2.5. This f ^ ^ ^ 

. • ,K .A = - 1) °' ^\ (24) 

IS m the expected range. \ „ ± r , 

^ I e^e^ , for t > T 

This shows that, after a time t ~ a few t, the effect of 

5.2 Role of external symmetry breaking perturbations 

the non-symmetric disturbance is exponentially amplified 

In reality, a strictly symmetric substrate, either spherical (e*/'^ > 1) by the instability mechanism, while the direct 

or cyhndrical, is impossible. Also, the chemical properties effect of the source is small in the sense that e^i < 1 

of the substrate surface are never perfectly homogeneous, even for t ~ few Tq. In this way, the symmetry break- 

A nominally spherical latex bead may contain a weak lo- ing mechanism manifests itself as an amplifier of small 

cal deviation of the surface curvature and a weak het- heterogeneous disturbance in the system, which can be 

erogeneity of the polymerization constant, kp, along the experimentally observable. 

surface. Thus we should suppose that there are a distur- An other way to think about the external perturbation 

bances which break externally the symmetry of the sys- is to define the time Tsq required for developing an of 

tern, and modify Ea. lfT!^ or its linearized form Ea. fltHl . a specified value e^: Equation leads to: 

A legitimate question is, therefore, "if and how the above / gS \ 

rs, = r,ln 1 + ^ . (25) 

symmetry breaking instability plays a role?" In short, the \ ^i / 

answer is that, despite these extrinsic factors, the insta- Changing the prescribed value or the external pertur- 

bility mechanism of symmetry breaking manifests itself in bation by orders of magnitude changes ts only by a 

the evolution of the gel's thickness, justifying our compar- small factor. This tells us that, as already announced the 

ison with experiments done in the above S iiS.ll We discuss scaling of the characteristic observable times is essentially 

it first in a formal manner, and then in the context of the given by Tg. It also tells us that the detailed knowledge of 

geometrical and chemical heterogeneities. the early dynamics is not essential in the definition of Tsq 

Within the linear approximation, the evolution equa- provided Tq is sufficiently larger than tq, so that a quasi 

tion Ea. H13|) is decomposed into the equation for each spherical state is obtained before the symmetry breaking 

mode, like Eas. (|15|) or H19|) . In the latter equation the process is observed. We know this to be true both from 

system's intrinsic heterogeneity may be represented as a our analysis and from experiment. 
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Now we describe how the parameter reflects the ef- standing the main feature of the dynamics if tq < Tq. We, 

feet of the heterogeneity of the surface curvature and of therefore, only mention about the restricted case where 

the polymerization rate. Under the linear approximation, those heterogeneities are switched on at a certain moment 

we only need to consider the profile of the substrate sur- of time after the symmetric steady state has been estab- 

face which can be described in terms of the radius ro(0o) s& lishcd. The result then reads 

a function of the angle e^: ro((?o) - ro + zi^*^°°^ cos(g0o). e^^^x / q 2 ^(gco) \ 

< = -7^ — — t ^^15— 9-1 • 27 

Additionally we consider the spatial distribution of the \ ^ / 

polymerization rate constant ~kp represented as a function The positive coefficient in front of A'^''^"'^ reflects the 

of 60: kpiOo) = kpo + A^t"""^ cos{q9o). A'f''°^ and z\(=''""' acceleration of the turnover of the gel material where ~kp is 

characterize the amplitudes of geometrical and chemical increased, while the minus sign in front of the second term 
perturbations, respectively. The geometric profile ro(0o) ^^"^ ^''^''^''^ Polymerization being slowed 

leads to the non-homogeneous curvature niOo), which has ^^^^^ ^^le^^e ^^"^ ^mf^'^e extrudes, or, where ^{60) > 
the following form, 

6 Discussion 



/v{gco) 

1 + ^ (q^ -l)cosiq9o) 

ro 



., ,11. r 1 1 - A 1- ■ Our analysis based on gel elasticity leads to the essential 
Along the hne of calculation m Appendix D, this expres- 

^, , 11, 1 ,r, -1- prediction that the spherical symmetry is always unsta- 
sion of the curvature should replace the factor ro in 

„ rm mi 1 ; ; 1 1 ; . c • tlc. Thc cxpcctcd scenario is that in a first step a quasi- 
-bjq.||34|). ihe normal stress on the substrate surface is 

, ^ . , spherical steady state is reached which should obey the 
therefore given by 

prediction contained in |S] and jl8|. Then on a time scale 

arrl =-k{0)T. (26) 

° significantly larger than the characteristic time for reach- 
As for the chemical heterogeneity in the polymerization ing the isotropic quasi-steady state, symmetry is broken, 
rate, kp{0o) should replace kp in Ea. l|13(l . In general these In the regime we discuss, governed by elasticity, these two 
effect can be summarized in the form of times are predicted to scale like the radius of the bead on 
— l3q'^'^'^"^\Aq'^'^'^"^^ /kp)+Pq^'^°\A'f'^°^ /tq), with dimcu- which the experiment is conducted. This scaling should be 
sionless numbers a,nd j3lf'^°^ . However, if these very robust in the elastic regime, since tq is the only length 
source terms have existed from the start of polymeriza- scale in the problem. In particular it should hold for wave- 
tion, the expressions of and j3^q'^°'^ are complex be- length larger than those considered here. All these expec- 
cause in the early stages of the gel growth none of the lin- tations are well born out by experiment 7 . In a number 
ear equation is valid. However, as we have already pointed of cases symmetry is not observed to be broken: this may 
out the exact knowledge of is not essential for under- be due to three different causes: - the experiment dura- 
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tion might not be long enough for the symmetry breaking the microtubule/centrosome (or microtubule/bead) sys- 
event to take place, - the gel/bead friction, considered in tem. The latter system does not exhibit an instability 
appendix E might further slow down the symmetry break- whereas the first does. The difference in behavior results 
ing process, - the gel might not behave fully elastically at from difference in boundary conditions. All these cases do 
very long time scales. In this latter case, a new time scale not consider the situation where filaments are crosslinked. 
would come into play, namely that over which a significant Actin gels are crosslinked and we propose that in two and 
stress may be maintained, and a new calculation should three dimensions these crosslinks change profoundly the 
be developed. We discuss various possible improvements behavior. Indeed, if the thickness of the gel layer is locally 
to our current analysis in appendix E. decreased, the compressive stress there, Crrl^Qi should ei- 

ther stay constant if full-slip boundary conditions are achieved, 
or decrease, irrespective of the thinning cause. The lateral 
displacement of the gel layer along the substrate surface 
might at most relax some of this local decrease of cTrrlrgi 
but it will never be able to increase it. Thus the mecha- 
nism described for non-crosslinked filaments do not 
apply to gels. 

A direct experimental assessment of the symmetry break- 
ing mechanism could involve monitoring simultaneously 
the depolymerization and the polymerization processes at 
the outer and inner gel surfaces, respectively. This is not 
an easy experiment. 



As explained in this manuscript the main ingredient for 
the occurrence of symmetry breaking comes from the ten- 
sile stress concentration where the gel thickness is small- 
est. This feature, added to a stress dependent depolymer- 
ization in the immediate vicinity of the gel outer surface, 
leads to an absolute instability of the system. This is in 
contrast with earlier models [SlllOj in which symmetry is 
broken at polymerizing gel bead surface. Their interpre- 
tation is most transparent in the one dimensional case; 
consider two opposing sides on which parallel filaments 
are grown. The force on individual filaments, i.e. the ra- 
tio of the total force (equal on both sides because of force 
balance) to the number of supporting filaments, is the key 
notion. The smaller the number of filaments participat- 
ing, the slower the effective polymerization rate; it is nat- 
ural to expect a force dependence and different scenarios 
have been discussed pill(J| if an unbalance between the two 
sides arises at some point, it grows since the "weak" side 
tends to become "weaker" . The two-dimensional version of 
this mechanism, simulated by van Oudenaarden et al. [U] 
is closely related to simulation and experiments done on 
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This pressure integrated over the surface, r^Af2 (3D) or 
rAd (2D), gives the total force exerted by this thin layer. 
These force in the respective dimensionality are shown on 
the right hand side of the equations below. Now, these 
forces must be the origin of the difference between the in- 
tegrated normal stress acting at r and that at r + Ar. The 
differences in the respective dimensionality are shown on 
the left hand side of the equations below. 

2(7 I I Ar 



Fig. 5. Forces acting on a curved slice of gel of thickness Ar 
and length rAO. The slice is under lateral tension because of 

forces a±±Ar. It is radially compressed because of the forces (r+Ar) Ail arr\r+Ar~i^ A[2 arr\r = f A[2 i(3D) 



{r + Ar)A9 and rAO ^ ^ ^ii^,(2D) 

r 

Appendix A. Heuristic derivation of the equa- Dividing the both hand sides of the above equations by 

.. £• u ■ I u I Ar, and letting Ar — > 0, we arrive at Eqs. IHl) and iffl. 

tions of mechanical balance & ^ h »_» 

Note that < for compressive stresses. 

Eqs. ||2Jl and (01 in §|3are the equation of mechanical bal- 
ance of stress components in the spherically and circularly 

Appendix B: Elementary physical mechanism 

symmetric geometries, expressed in respective relevant co- 

of the stress concentration 

ordinate systems. Instead of deriving these from the famil- 
iar form in the Cartesian coordinates (symbolically writ- We describe the basic mechanism of the stress concentra- 
ten as V • CT = 0) through coordinate transformations, we tion by an illustrating example with a very simple geom- 
will present an elementary physical interpretation of the etry (sec, Fig|6Ka)). 

equations of mechanical balance. It might help to under- Suppose that there is a long elastic rod whose diame- 

stand how the lateral tensile stress and the normal com- ter d{z) is inhomogeneous along its long axis, z. We now 

pressive stress are related. See, FiglSj Consider, within apply a tensile force to this rod by pulling its ends apart, 

a layer of actin gel occupying the radii tq and + h, a Once the balance of force is reestablished within the rod, 

slice of gel between the radii r and r + Ar spanning a the total tensile force integrated over a sectional plane 

solid angle Af2 (3D) or an angle 6 (2D). The lateral ten- perpendicular to the z-axis is constant along of the z- 

sion a±± gives an effective surface tension AF = a±j_Ar coordinate. Thus the tensile stress azz averaged over this 

to this slice. Because of the curvature radius, r, of this section is inversely proportional to its area, 7r(d(z)/2)2. By 

slice, a sort of the Laplace pressure, 2Ar/r (3D) or AF/r such geometrical effect, the tensile stress is concentrated 

(2D) is generated towards the center (r = 0) of curvature, at the thinnest part of the rod. 
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(ex. the momentum, the charge, etc.) which is confined 
along some direction(s). 



Appendix C 



Fig. 6. Stress distribution in elastic structures of varying 
thickness; (a) Elastic rod of thickness d{z) under tension due 
to forces F acting at the ends. The tension at coordinate z: 
a^^ ~ F/A{d), where A{d) = jd{zf. (b) Analo gous situation 
in a gel layer of thickness h{Oo) under integrated tension T. 
The tensile stress can be approximated as ~ T^J)- 

We could mention a analogous situation in an electric 
wire transporting a steady electric current. If the thick- 
ness of the wire is inhomogeneous, the electronic current 
density is high in the region where the wire is thin, by the 
same geometrical effect. The electric current density plays 
the role of the tensile stress (Jzz in the former case. In fact, 
the stress is the current density of the momentum |12| . 

We may compare these quasi one dimensional exam- 
ples with the geometry studied in ^ (see Fig|Sfb)). In 
the latter situation, the gel is under lateral tensile force T. 
By the same reasoning as above, the lateral tensile stress 
is large where the thickness is small. Although the shear 
force between the gel and the cylinder would weaken this 
effect, the basic mechanism still works. 

The geometrical effect discussed here is quite universal: 
We only need a current density of some physical quantity 



We demonstrate that the perturbations of the surface pro- 
file have a limited influence on the stress, practically con- 
fined within a "skin depth" near the outer surface of the 
gel layer, where the skin depth is of the order of the wave- 
length of the perturbation (see Fig|31in the text). 

(i) To the spherically symmetric gel layer, we introduce a 
xyz-coordinate system so that its xy-coordinate plane is 
tangent to the outer surface of the gel layer at its origin, 
X — y = z — 0. We define the sign of ^-coordinate so that 
the bulk of the gel is on the side of z < 0. We will consider 
a small neighborhood of the origin so that the curvature of 
the gel surface is negligible. This apparently flat gel layer 
is under lateral tension along the xy-plane. 

(ii) We introduce a slight sinusoidal perturbation of the 
surface proflle of the gel layer, without allowing the dis- 
placement of the gel material. The perturbed surface pro- 
file is written as 

z = eg i?e[e*'5^+'*], with the amplitude eg, the wave num- 
ber Q and the phase (j) being constant. 
(Hi) We then let the gel to relax until the mechanical 
balance is reestablished within the layer. By this process 
the stress components aaf3 with a ^ x, y, or z are also 
perturbed. We denote by Saafs the perturbed part of the 
stress components. These Saaf3 obeys the following equa- 
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tions: with the form e^^ for z < 0. This indicates that the influ- 

d d d 

-rr-Suxa + -K~5'^va + TT^'^za — Oj encB of the perturbations to the surface profile with the 



wavelength ~ is practically limited within a region 
with a "skin depth" ^ Q^^ from the outer surface of the 
gel layer. 



OX ay oz 

with a = X, y, or z We assume that the usual linear elas- 
ticity relationship applies to the system. Then the per 
turbed stress components is related with the displace- 
ments {ux^ Uy,Uz) from the unperturbed state through the 

equation: Appendix D. The derivation of Eqs. O llj t and 

7 7 

-with Here we show how the stress distribution within the gel 

u a = Ij'^"" -I- ^HH) layer is calculated for the model described by (a)-(c) in 

fj^l In the polar coordinate, the equations of mechanical 

where /i and K are the shear and bulk moduli, the surhces 

„ , , , r 1 r -1 mi balance in the gel layer is written as follows: 

a, p and 7 take x, y or z, and {Xx, Xy, Xz] = (x, y, z\. ihe 

summation index 7 runs over x, y and z, and 5^/3 is the d ^ \ , ^ 

Kronecker's delta. 

(iv) The question is how the quantities Suap depend on z 

for z < 0. In the lowest order of eg, Saap, and therefore We integrate the left hand side of these equations with 

the displacements u„ should depend on x sinusoidally with respect to r from vq to ro + h, noting the boundary con- 

the wavenumber, Q. Therefore, the above equations can ditions Eqs.UHI) and ©. The result reads 
be reduced to the following matrix equation: ^ 

T - —f = -roCTrrlm 



Q^{rarr) + T^cTre- (Tee = 0, 
d d 

^(w^e) + ^(Tee + (Tre = 0. (28) 




\^(K+|M)ig^ {K+fi,)^-2t,Q^j \u,j \0j and 

- d 

This equation can be finally reduced to the following T + tttT = 0, 

Oo 

equations: , i • 1 1 i 1 • m rra+h , 

where we have introduced the total tension, J = J^^^ crggdr^ 

{dz Q ) "0 analogue for the shear stress, T = J^°^'^ cTre dr. 

where V' is a combination of and Uz- From this equa- Based on our estimates of the shear stress, (7^9 '~ pi-^qih* / rof' ^ 

tion, we find that the displacements should depend ex- and the perturbed part of the tensile stress, 5agg ~ eqBh*/ro 

ponentially on z. Among mathematically possible forms (see we can evaluate the terms on the left hand side of 

g±Q2^ we discard the form e"*^^ since this factor grows the above equations. Since ^ ~ 1 under our limitation of 

exponentially towards negative z axis. We are then left the wavenumber q, we have, ST ~ - Beqh*'^/ro and 
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f ~ jleqh*^ /ro^ , where 5T is the perturbed part of 

the total tension T . Assuming jl ^ B, we find that T and 
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09 



T are smaller than T and by a factor of h* /tq- We, 



therefore, ignore the terms with T, and have the following 
equations: 



(29) 



The first equation requires the lateral balance of the inte- 
grated tension T, while the second equation requires the 
homogeneity of the normal compressive stress on the sub- 



strate surface, arr\ro- 



From (|29|l we can calculate aee 



ro+h 



, and — 



We employ the "stacked rubber band model" |51llllll2| 

for the lateral tensile stress aee, as we did for the sym- 
metric case (see Eq.([Pl in the text). Here we take into 
account the possible lateral displacement of the gel layer 
upon the reestablishmcnt of the mechanical balance. We 
introduce an unknown function 9{9o) such that the ma- 
terial of gel layer originally at Oa is moved to 9 (da) upon 
the reestablishmcnt of the mechanical balance (see Figl7|) • 
The elongation ratio, (r— ro)/7-o in the Eq.(^ is, therefore, 
replaced by the form which dependents on the parame- 
ter Oq: {rd6{9o) ~ rod9o)/rod9Q. Thus the lateral tension 
is written as 



(Tee = B 



ro d9o 



- 1 



(30) 



where, the shear deformation within the layer has been 
consistently ignored. (The justification of this approxima- 
tion concerning another source of error will be discussed 
in the Discussion section, see Appendix E.)With the dcfi- 



Fig. 7. (a) 9o is defined as the angle with respect to a reference 
line, when a material point (for example, the black dot) is 
located before the gel is deformed, (b) as a result of elastic 
deformation, the material point characterized by 6*0 is displaced 
to a new position at 9. The function 9{9o) characterizes the 
elastic deformation. 



nition of T given above, we obtain 



T = B 



hi9o) 



h(9ol\ d9{9^ 
2ro J d9o 



- h{9o) 



(31) 



The function 9{9o) can be related to h{9o) through the 
first equation in l|29|) . which requires that T is constant. 
To fix the value of the constant, T, we recall an apparent 



condition J^- 



27r de(eo) M - 



de, 



d9o = 271. The resuh of T is 



T = B 



27r 



d9o 



h(eoY 



^^(^0)+ 2ro J 

2tt- 



hi9o) 



H9o) + !^ 



-d9o 



(32) 



den 



can thus be finally determined in terms of /i(0o) 



(which we do not show explicitly) . From the second equa- 
tion of and from Ea.l|5n|l. we have 



cree\ra+h{eo) 

T I 1 



H9o) 
■ h{9, 



hi9or] 


-1 









2ro 



T 



(33) 
(34) 
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To reach the expressions Eas. (|llll and p2() in the text, we bility result. For the other modes, especially for q < ^f?-, 
may simply substitute the form /i (00 ) — h* [1 + cos((7 0o)] i we should take into account the friction on this surface 
into Eas. (|32|) - (|33|l . and develop them with respect to up due to the temporal linkage between the actin filament 
to the linear order. with the substrate ([ill). (About the discussion of the re- 

lation between the friction and the temporal linkage, see 
Appendix E |19II15| .1 As a modification of the present model, we have 

incorporated the finite friction force on the substrate sur- 



Below we mention briefly aspects which could be improved 
in our present analysis: 

1. Extend the analysis to modes with q ^ ^^r^' 
Although the instability against the disturbances of modes 
q ~ is sufficient to destroy the system's stability, our 
analysis can say nothing about what is the fastest, or the 
most unstable, mode of the disturbance. The fact that 



face, which is proportional to the slipping velocity of the 
gel along the surface. Though details will not be shown 
|13| . the result indicate that, while the all modes remain 
unstable, the instability is weakened for long wavelengths, 
i.e. for small values of q. 

3. Extend the analysis to non-linear regime: 

^, , i • i- i- r ii i 1 1 • • Our analysis does not infer how the comet of actin gel is 

the characteristic time of the symmetry breaking m our ^ 

, . 111 X ii i ii ii formed and continues to grow after the symmetric shape of 

analyses gives reasonable values suggests that the other ^ j t- 

J r i u • 1 i T J i J- the layer around the bead is lost. This is a nonlinear prob- 

modes of perturbations might grow, if they do, at a rate 

.i-i, j-ii 1 ijTTi lem. ini have demonstrated in their numerical modelling 

not highly exceeding the one we have analyzed, m fact, ° 

a- ^ , /2^i i 1- / iiN- that the comet formation shows its optimal performances 

some cftorts to renne the present analysis (see below) m- ^ ^ 

J. , , r , 27rrn ii ■ i 1 -i-i ■ 1 J for a certain parameter value related to the depolymeriza- 
dicate that, for q < =tt^, the instability is weakened or 



h' 

27rrn 



1 1 .1 r ii 1 -ZTTm ii • tion at the substrate surface. The comparative study from 

even suppressed, while for the modes, q > ^j^r-, there is t- j 

r • ui J J TT iu J our point of view is yet to be done. (As for the steady 

no sign of appreciable q dependences. However, the mod- r- j \ j 

1 ^. r ii • 1 1 1 i ii 1 • r growth of the comet from Listeria, see (HTl)). 

ulation of the micro-scale comparable to the mesh size of ^ \t—i" 

the gel is not accessible by the continuum approach. (Note 4- Extend the analysis to soft beads: 

that the quasi one-dimensional analysis of | 1(J| claims the Endosomes, lysosomes, vesicles and fluid drops deform as 

instability of the mode q — 2.) the comet develops, revealing the importance of mechan- 

2. Remove the full slip boundary condition: ical stresses pil2()ll21| . The deformation of a fluid drop 

As discussed in S]3| we have justified this boundary con- has been fully analyzed within the framework of the elas- 

dition when we analyzed the evolution of the modes q ~ tic analysis and shown to be quantitatively in agreement 

2^?- , since, there, the choice of boundary condition on the with the experiment |22j. The symmetry breaking onset 

substrate surface is expected to be insensitive to the sta- remains to be worked out. As proposed in the discussion. 
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the early stages would discriminate between the different 
mechanisms. 

5. Removal of the assumption of the isotropic gel: 
Generally speaking, the micro-structure of the gel poly- 
merized from a surface must distinguish the radial direc- 
tion from the lateral ones. Especially, for the actin gel 

branched by the help of the protein Arp2/3 is shown to „. „ „ „ , . i r . ,. , , 

' 1? ig. 8. Generation oi a mechanical trustration caused by 

have a topology like a "forest" rather than the network gi^juitaneous polymerization and lateral deformation a rect- 

[?]. It will be the entanglement among the branches of .^^^^^^ ^-^^^ (marked block) is displaced and deformed 

the "trees" of semi-fiexible filaments that supports the while keeping the connectivity with its neighboring piece (block 

tensile stress within the gel. Though we expect no quali- shown by dotted lines), 
tative change of our result upon the incorporation of the 

elastic anisotropy of the gel, there should be quantitative tion is advanced in a portion of gel layer than elsewhere, 

differences. For the further analysis, we also need the ex- the stress concentration is most likely to occur and the 

perimcntal data on the anisotropic elastic constants T7. degradation will further accelerated there. Visco-elastic or 

„ „ , , ,, , . , ,, , , ., . ,. „ frictional effects in the bulk gel or on the substrate sur- 

6. hxtend the analysis where the gel density is spatially 

, , face, respectively, may limit this positive feedback loop. 
heterogeneous: 

mi £c i rxi i- 1 1 i -i r ii j. t j.- The detailed discussiou wiU bc the task of futurc works. 
1 he eftect of the spatial heterogeneity of the catalytic ac- 
tivity of enzyme may have several aspects. The one which 7. Take account of the frustration of stress in gel: 
has been discussed in ^ ^5.2\ is the the modulation of the In the analysis of f^l the total tension T = J^°^'^ aggdr 
polymerization rate, kp. The other aspect which is related has been calculated by substituting the expression of the 
to the spatial heterogeneity of elastic moduli of gel may stacked rubber band model, Ea. (|30|l . This operation ig- 
also deserve consideration. In fact, the heterogeneity of nores the fact that the gel material at two different radii 
the elastic moduli will be closely related to the hetero- are created at different points of time. The Fig|Slillustrates 
geneity of the factors k^ and Cd both concerning the de- how the simultaneous polymerization and lateral deforma- 
polymerization processes kinematically and energetically, tion create a mechanical frustration within the gel mate- 
respectively. It is therefore impossible to predict where rial: (a) Consider a thin slice of actin gel created at the 
does the thinning of the gel layer proceed most rapidly, substrate surface during a short time interval, say, between 
However, the rule of thumbs is again that the positive ti and ti + dt (the dark gray region occupying the angle 
feedback mechanism mentioned above: once the degrada- ^do). We may expect that this part of gel which is just 
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grown bears no lateral stress, to a good approximation. In f = 30nm for the mesh size of the actin gel. We, there- 

the context concerning this slice just above the substrate fore, suppose that our approximation is pretty good for 

surface, we would then set dd{9o)/d9Q = 1 in Ea. H3U|) . (b) the above parameter range. 

After the consecutive time interval, ti + dt < t < ti + 2dt, 

the same spatial region, which is now indicated by the 

dotted lines, is occupied by a newly grown gel under no 

lateral tension. Thus again d9{dQ)/dda — 1 for this region. 

However, as for the previously grown material which we ^- Verbal, F., V. Laurent, A. Ott, P. Chaikm, J. Prost. Mea- 
have marked in dark gray in (a) , it now occupies the region surement of the elasticity of the actin tail of Usterta mono- 

. , . J ^, • • 1 /I • ■ J 1 \ cj/togenes Eur. Biophys. J. 29 (2000) 134-140. 

just outside the origmai one (shown again m dark gray), 

and occupies the angle AO. Generally, AO is different from 2- Taunton, J., B.A. Rowning, M.L. Coughlin, M. Wu, 

An r ,1 • iri R.T. Moon, T.J. Mitchison, C.A. Larabell. Actin-dependent 

Avo as tar as there is a global lateral displacement oi gel ' ' ^ 

, . , 1 , . . , , , , , , , r, , , mi • 1 Propulsion of Endosomes and Lysosomes by Recruitment of 

during the time interval ti + dt < t < ti + 2dt. 1 hus m the 

N-WASP. J. Cell Biol. 148 (2000) 519-530 

context concerning this region in dark gray, we would set 

^ I -in / 1 • rm, mi • , 1- ■• • 1- . 3. Loisel, T.P., R. Bouiemaa, D. Pantaloni, and M.F. Car- 
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